d01akf

d01akf © Numerical Algorithms Group, 2002.

Purpose

D01AKF 1-D quadrature, adaptive, finite interval, method suitable for oscillating functions

Synopsis

[result,abserr,w,iw,ifail] = d01akf(f,a,b,epsabs,epsrel<,lw,liw,ifail>)

Description

 
 D01AKF is an adaptive routine, using the Gauss 30-point and 
 Kronrod 61-point rules. 

 Because this routine is based on integration rules of high order,
 it is especially suitable for non-singular oscillating 
 integrands.
 
 D01AKF requires the user to supply a function to evaluate the 
 integrand at a single point.

Parameters

d01akf

Required Input Arguments:

f                                     function (User-Supplied)
a                                     real
b                                     real
epsabs                                real
epsrel                                real

Optional Input Arguments:                       <Default>

lw                                    integer  d01akf09
liw                                   integer  d01akf11(lw)
ifail                                 integer  -1

Output Arguments:

result                                real
abserr                                real
w (lw)                                real
iw (liw)                              integer
ifail                                 integer